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ABSTRACT 

We present a method to extract the redshift-space distortions (3 parameter in configuration space 
with a minimal set of cosmological assumptions. We show that a novel combination of the observed 
monopole and quadrupole correlation functions can remove efficiently the impact of mild non linearities 
and redshift errors. The method offers a series of convenient properties: it does not depend on 
the theoretical linear correlation function, the mean galaxy density is irrelevant, only convolutions 
are used, there is no explicit dependence on linear bias. Analyses based on dark matter iV-body 
simulations and Fisher matrix demonstrate that errors of a few percent on /3 are possible with a full 
sky, 1 [hr^ Gpc)"^ survey centered at a redshift of unity and with negligible shot noise. We also find 
a baryonic feature in the normalized quadrupole in configuration space that should complicate the 
extraction of the growth parameter from the linear theory asymptote, but that does not have a major 
impact with our method. 

Subject headings: cosmological parameters - cosmology: observations - cosmology: theory - dark 
energy - large-scale structure of universe - methods: statistical 



1. INTRODUCTION 

The gravitational instability picture involves the 
growth of early time primordial fluctuations into the ac- 
tual large-scale structure observed through the galaxy 
distribution. This growth depends in principle on the 
underlying theory of gravity and the cosmic expansion 
history. It is therefore important to be able to measure 
the growth history to obtain useful cosmological infor- 
mations. 

One way of determining the growth of structure is 
through the apparent anisotropy of the galaxy distribu- 
tion in redshift space, caused by the line-of-sight (LOS) 
component of the galaxies peculiar velocities. On large 
scales, under the linear perturbations regime, the two- 
point correlation function presents a squashing along the 
LOS and correspondingly the power spectrum appears 
enhanced for wavevectors directed along the LOS ( Kaiser 



19871). The anisotropics are governed by the parameter 
jj that depends on the growth function and galaxy bias. 

Redshift-space distortions (RSDs) h ave been the sub - 
ject of many analyses, as reviewed in Hamilton (1998). 
Examples of recent studies involving the latest large sur- 
veys are the following. The Two-Degree Field Galaxy 
Redshift Survey enabled the RSD measurements on the 
correlation function (Peacock ct al. 2001; Hawkins et al. 
[2003J and power spectrum ( jPcrcival et al. 2004j . I'he 
Sloan Digital Sky Survey permitted ot her RSD mea- 
surements on the correlation functi on (iZehavi et al.l 
2005||Okumura et al | 2 008 | |Cabre fcG aztaiiaga 2009a b) 
and power spectrum ( 'I'egmark et ar||2004[ |2006| ). The 
VIMO S-VLT Deep Survey and the 'l^YJm Surrey were 
used in Guzzo et al. (2008) for RSD determinations from 
the correlation function. After this work was submit- 
ted RSD s tudies on the WiggleZ and BOSS ca talogs also 
appeared ( [Blake et al.||201l[ |Reid et al.||2012[ ). 

Since the linear theory description starts to be valid 
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only at very large scales, and lacking a complete model 
for general non linear cosmic fluctuations, an extension of 
the theoretical description has been attempted to non lin- 
ear and quasi-linear scales thanks to em pirical method s 
based on the so-called streaming model (Peebles 1980), 
consisting of linear theory and a convolution on the LOS 
with a velocity distribution. The model was first adopted 
on small scales and highly non linear regime to describe 
the fingers-of-God (FOG) elongation alo ng the LOS du e 
to random motions of virialized objects ( Jackson]|1972 ). 
Also, fitting functions ba sed on simulation resu l ts have 
been used, f or ins tance, i n Hatton fc Cole ( 1999), Tinker 
eral]([2006| and [Tin ker' ("20071). Recently, iT was shown 
Dy Percival & White ( 2009 ) that on quasi- linear scales a 



streaming model with a Gaussian velocity dispersion is 
a good general fit to the redshift-space power spectrum 
The goodness of the streami ng model wa s also demon- 
strated lately, for example, in'Guzzo et al.' ('20081), "Cabre 
& Gaztaiiaga (2009a), Blake ct al. (2011), Beutlcr ct al 
l213l2| ), |Marufii et al.| ( |2012^ ) and |Chuang fc Wang, ( |2011f 
We caution that the model breaks down at small scales. 
This happens even for an unbiased dark matter model 
(see, for example, Taruya et al. (2010)) and complex 



galaxy bias issues niaybe required to be u n der strict con 



trol (li ke shown in Okumura & Jing ( 2011 ); Reid & White 
(2011J). Also, how modes representing large scales are 
1( 



sam pled in simulations finite volume boxes might play a 
role ( [Matsubara 2008} [Taruya et al.[[2009[ ), especially if 
their eflect is ot the order of a few percent and in the 
presence of bias, that might further reduce sampling ef- 
ficiency. Nonetheless, being aware that there exist in the 
literature works that claim either a good or not so good 
performance of the streaming model, and that anyway 
often use it as benchmark, due to its theoretical simplic- 
ity and appeal, in this work we will rely on a generalized 
version of it, with the crucial limitation to deal with ap- 
propriately large scales. 
In this paper we find a novel way to measure the (3 



2 



parameter on quasi-linear scales in configuration space 
without utilizing a theoretical description of the linear 
correlation functions. We generalize a sim ilar result in 
Fourier space obtained in Cole et al. ( 1995 1 . Our method 
does not depend on the galaxies mean density and does 
not employ deconvolutions but only convolutions, an ad- 
vantage in terms of stability. 

We also note the presence of a baryonic feature at 
around 100 Mpc in the so-called normalized quadrupole 
in configuration space. We arc not awarc[^ of other pa- 
pers menti oning it , apart from a hint from Fourier space 
provided in |Taruya et al.|(|26ld | ) and in relat ion to biasing 
from density peaks in Desjacques & Sheth (2010). Prob- 
ably, the reasons for this neglect are that tne pioneering 
RSD papers did not contain a baryonic component and 
subsequent theoretical analyses were concentrated on rel- 
atively smaller scales that were better sampled by the 
surveys being studied. 

The plan of our work is as follows. We start in Sec- 
tion [2] with a review of the theory of RSD in Fourier and 
configuration space in the linear regime and we discuss 
the streaming model extension. In Section [3] we describe 
our method and the N-body simulations that we use to 
test it. Section |4] demonstrates how stringently the pa- 
rameter j3 can be measured, thanks to a full statistical 
analysis with the simulations and a Fisher matrix study. 
We conclude in Section [Sl 

2. THEORETICAL BACKGROUND 

In what follows we provide a theoretical background 
both in Fourier and configuration space. It is convenient 
to be able to switch between the two points of view be- 
cause, although they are theoretically equivalent, when 
dealing with real data, they are confronted in practice 
by somewhat different systematics. Fourier space has an 
advantage regarding the theoretical treatment because in 
linear theory the Fourier modes are independent. How- 
ever, it has been noted that perturbation t heory seems to 



be more e ffective in configuration space (Sanchez et al. 
[2008 ; Matsubara|2008[|Taruyaet al |2009[ ), i.e. that con- 



vergence happens to be taster with respect to Fourier 
space. 

We want to stress that in this work the Fourier treat- 
ment is utilized only to derive convenient expressions 
in configuration space. The observational inputs of our 
method and manipulations upon them are just in config- 
uration space. 

2.1. Theory in Fourier Space 



It was shown by Kaiser (1987) that in the large-scale 
linear regime and in the distant-observer approximation 
the coherent gravitational infall modifies the redshift- 
space power spectrum in Fourier space: 



(1) 



where P{k) is the real space power spectrum of the parti- 
cles that trace the density fluctuations 5, Hk is the cosine 
of the angle between the wavevector and the LOS, the 
subscript s indicates redshift space, and /3 is related to 



^ After this work was submitted, [Samushia et al.] | |2012| showed 
a plot where the feature is visible in simulations results. Also, Eyal 
Kazin(2011) in a private communication sent us and commented 
on a similar plot. 



the perturbations growth rate in linear theory. We define 
the linear bias b as the ratio of the tracer overdensity S 
with the matter density 6m, 



(2) 



The relation between the matter and tracer power spec- 
trum is simply 

P{k) = b'Pmik). (3) 

We will assume for simplicity a constant bias since we 
are dealing with large scales and will discuss later the 
possibility of eventual departures. The crucial parameter 
P is 

/(a„) ^I dXnP 

b b d\na ^ ' 

in which a is the cosmic scale factor and D is the lin- 
ear density growth factor. The variable / can be well 
approximated by 



(5) 



with 7 :^ 0.55 if general relativity is assumed to be valid 
( Linder|2005 ); flmia) is the fractional density of matter: 



m{a) 



(6) 



where fiom is the matter density today (the subindex 
always indicates actual quantities) and 



H{a) = Ho^J^oma-^ + ^ka'^ + r!odca-3(i+-) (7) 

is the Friedman equation that defines the Hubble param- 
eter H{a) in which figde and w are dark energy density 
parameter and equation of state, and fi^ is the curvature 
energy density. 

The linear power spectrum can be decomposed in a 
multipole expansion ( [Cole et aL]|1994[ ) 

P,(fc, ^ik) - Po{k)ro{Hk) + J^2(s)p2(/ifc) + n(s)n(Mfc), 

in which the multipole coefficients are expressed as 



21 + 1 
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d^kVii^J-k) Ps{k,fJ.k), (9) 



where Vi{^ik) is a Legendre polynomial. Each of the co- 
efficients is given in terms of the real space tracer power 
spectrum 



Fo(fc) = ( 1 
P2{k) 
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P(fc), 



Pi{k) 



8^ 
35 



P{k). 



(10) 



(11) 



(12) 



Now, if the ratio of the quadrupole to the monopole is 
taken, this is expressed by a rational function of the /? 
parameter ( Cole et aL]|1994 ): 



Q(fc) = 



Po{k) 



1 + u + 



(13) 
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Note that, as far as linear theory is vahd, this ratio, that 
we call normalized quadrupole, is scale independent. 

2.2. Theory in Configuration Space 

It is possible to transfer Kaiser's treatm ent from 
Fourier to configuration space ( Hamilto"n||1992 ), 

tt) = Us)Voi^^) + 6(s)P2(m) + Us)V^{^i), (14) 

in which tt is the distance separation along the LOS and 
a is the perpendicular separation. The absolute distance 
of separation is indicated by s = Vc^ + tt^, ^ = n/s is 
the cosine of the angle between the separation vector and 
the LOS. ^ 

The multipoles of ^(tr, tt) rl are given in terms of the 
real space correlation function £,{r): 



Us) = { 1 

6(5) 
Us) = 



3 
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[U)~U)], 
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and 



U) 



U) + 2^(0 - 



^{r)r^dr\ 



(15) 
(16) 
(17) 

(18) 
(19) 



If the real space correlation function ^(r) is computed 
from linear theory, these equations describe the typical 
squashing of the correlation function along the LOS of 
the Kaiser effect. The multipoles of £,{(j,7t) can be ex- 
tracted through the projections 



U) = ^ / UY^dr' 



Us) 



21 + 1 



-1 



Hamilton ( 1992 1 



1-1 
showed 



acT,n)Vi{fi)dii. (20) 
that the normalized 



quadrupole, defined as 

Q{s) 



Us) 



i3/s^)J^Us')s'^ds'-^o{sy 



(21) 



can be used to estimate the /3 parameter, since the Kaiser 
result, valid in linear theory and large scales, brings to 
the constant ratio 



Qis) 



1/3- 



l + |/3 + i/32- 



(22) 



All these equations mimic the Fourier description, in 
particular in both cases the normalized quadrupole on 
the large scales described by linear theory is a constant 
plateau. 

The multipole correlation functions are related to the 
corresponding power spectra by 



Us) 



■I 



1 
2"^ 



dke ii{ks)Pi{k), 



(23) 







^ |Kazin et al.| l |2012| have shown an alternative clustering wedges 
decomposition. 



here ji indicates the spherical Bessel functio n of the Ith. 
order and the following identities presented in Cole et al. 
( 1994 1 provide useful connections, that will prove helpful 
also later on. 



Us)-Us)^- 



1 
2^ 



dk j2iks) Pa{k) (24) 



Us) + lUs)-lUs) 



1 

2tt^ 



dk k^ ji{ks)PQ{k). 

(25) 



2.3. Departures from Linear Theory 

On very small scales the random motions of the galax- 
ies will produce the FOG phenomenon: an elongation 
of the correlation function mostly on the LOS direction. 
The FOG effect can be mimicked by a convolution of the 
correlation function ^(cr, tt) with the distributi on func- 
tion e ither of random pairwise velocities to give ( Peebles 
T980|: 



C(f^,7r)=/ i{a,7T -v)f{v)dv 



(26) 



where the peculiar velocities are expressed in comoving 
coordinates. The convolution device is usually defined as 
the streaming model. 

The random motions of particles belonging to a virial- 
ized system depend on how the galaxies are selected and 
have been represented in the literature by an exponential 
form, 



f{v) = 
or a Gaussian form 

/(«)- 



1 



(Ty^/2 
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V2\v\ 



V 

2^ 



(27) 



(28) 



where is the pairwise peculiar velocity dispersion. The 
exponential and Gaussian forms have proved to be good 
fits of the observed data and simulations (see , for ex- 
ample, Peacock fc Do dds 1996; Hamilton 1998} Peacock 
et al.|2001||Guzzo et al |2 008; Cabrc fc Gaztana ga"2009a| 



Blake et a 



120111 I Beutler et al. 2012,|Marum et al. 2012 



Chuang & Wang 2 01 ip . In this work our main concern 



IS on relatively large scales and our results do not de 
pend much on the choice of the distribution function. 
An important addition to the applicability of the abov e 



convolution has been offered by Percival & White (|2009|) 
They showed that the convolution with ay as a free pa- 
rameter could be useful to describe quasi-linear scales 
important for the transition to linear scales. As we men- 
tioned in the Introduction, at small enough scales non 
linearities and galaxy bias will start to spoil the simplic- 
ity of the model. How small these limit scales have to 
be depends also on the bias of the set of galaxies under 
study. 

Carrying over the convolution to Fourier space, one 
simply gets an additional factor to the Kaiser formula- 
tion, and what we define as the observed power spectrum 
looks like 



P°^^(fc,/Xfc) = {l + l34)'Pik)f{aM, 



(29) 
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where f{<7ykfik) is the Fourier transform of the velocity 
distribution. If the distribution is a Gaussian, it foUows 
that 

/KfcMfe) = e-(^"'=^'=)'/2, (30) 



while if it is an exponential we have that 



(31) 



Both forms have the same Taylor expansion for a small 
(jy up to second order. The projected multipoles of the 
observed power spectrum, that we also refer to as the 
observed multipoles, can be found from 



prik) 



21 + 1 



dtikVMP°'''ik,tik)- (32) 



Eventual random errors on redshift determinations can 
be described by a convolution with a similar Gaussian. 
Given the literature results favorable to the streaming 
model, we assume that as long as we refer to scales larger 
than the effective parameter, all the above mentioned 
deviations from linear theory with a linear bias can be 
reasonably controlled by the streaming model with the 
free parameter a^. Clearly, if one is able to predict the 
dispersion parameter from first principles, all the atten- 
tion could then be concentrated on the /3 parameter, that 
could be determined with more precision. 

3. A NEW METHOD TO ESTIMATE THE 
REDSHIFT-SPACE DISTORTIONS PARAMETER 

We now introduce our new method starting from the 
problems that origin ate from the normalized quadrupole 
defined in Equation ( 21 ). We test our method with the 
resu l ts fro m the A^-bocIy simulations by |Sato fc Matsub^ 



ara 



(2011) to estimate the level of systematic and sta- 
tistical errors for the recovery of the parameter /3. We 
also perform a Fisher matrix analysis to check on the 
expected error magnitude. 

As in the Sato and Matsubara simulations, our chosen 
cosmological parameters are: f^A = 0.735, = 1 ~ ^Ai 
= 0.0448, w = -I, h = 0.71, Us = 0.963, and 
CTg = 0.8. These are, respectively, the fractional density 
of dark energy, dark matter and baryons, the dark en- 
ergy equation of state, the Hubble parameter (in units of 
lOOkms""'^ Mpc~^), the primordial perturbations spec- 
tral index and the normalization factor relative to the 
variance of the density fiuctuations contained within a 
scale of 8/i~^Mpc. The 30 simulations that we utilize 



were done with the Gadget2 (Springel 2005) code with 
Np — 1024^^ dark matter particles, a box size Lbox ~ 
1000 /i^^ Mpc, a softening length = 50/i~^kpc, an 
initial redshift Zjni = 36, an output redshift z = 1. 

As mentioned by the authors there might be a slight 
systematic effect due to finite-mode sampling in a limited 
number of finte-size boxes. To quantify such an uncer- 
tainty we set up for an approximate correction t o the 
data. As sug g ested in Sato & Matsubara (2011) and 



Taruya et al.| (|2009| , we correct for finite- mode sam- 
pling by using simulations at a redshift when linear the- 
ory dominates. The simulations at the initial redshift 
Zini = 36, also kindly provided to us by Sato and Mat- 
subara, are ideal for our scope. We proceed by rescaling 



the Zini — 36 monopoles and quadrupoles to z = 1 us- 
ing linear theory. We then subtract those extrapolated 
monopoles and quadrupoles from the original z = 1 snap- 
shots. We thus obtain at this point a representation of 
only the non linear contributions due to the fact that the 
rescaled functions have the same initial random seeds as 
the original z — 1 snapshots. As a final step, we gener- 
ate random linear monopoles and quadrupoles, properly 
correlated among each other, and we add them to the 
above non linear-only parts. We have checked, using the 
untouched z = I simulations, that the final statistical er- 
rors on /3 do not change significantly, therefore justifying 
our approximate correction. We find that the extraction 
of the /3 parameter is corrected from a bias of about 4%- 
5% with respect to carrying out a parallel analysis on the 
uncorrected data. This might be due to the fact that the 
quadrupole, that is very sensitive to the /3 parameter, is 
quite susceptible to the finite-mode sampling problem, 
since it is, for example, less sampled than the monopole. 
Since we are expecting statistical errors on /3 at the few 
percent level, as the simulations and Fisher matrix anal- 
yses exposed later show, we employ our corrections in 
the rest of the paper. 

In Figure (fl]) we show the normalized quadrupole de- 
rived usi ng the linear pow er spectrum computed by 
CAMB (Lewis et al. 2000) for the fiducial cosmology, 
that implies tor dark matter /3 = 0.85, and a Gaussian 
velocity dispersion with a^^ — 5.7 h^^ Mpc. It can be no- 
ticed that the asymptotic constant form is reached only 
at very large distances, where the error due to cosmic 
variance is supposed to be important, and residual bary- 
onic acoustic oscillations at scales around 100 Mpc 
can potentially con fuse the estimation of the plateau. 



Taruya et al. (20101 find oscillations in the Fourier space 
ratio of the quadrupole and monopole power spectra. We 
speculate that the baryonic spike we observed using a 
simple phenomenological streaming model will b e some 
what enhanc ed due to the extra terms found in Taruya 
et al. (2010) for unbiased dark matter. The feature 
should perturb the extraction of the growth parameter 
in configuration space from the asymptotic limit of the 
normalized quadrupole if scales close to its location are 
employed. 

One could proceed by comparing the observed nor- 
malized quadrupole with the theoretical estimated one. 
However this would be done at the cost of introducing a 
dependence on the cosmological parameters, which would 
have to be marginalized over to obtain a statistically in- 
dependent error on /3. The above marginalizations are 
going to inflate the error on /3. If not predicted properly, 
the dispersion has to be marginalized over. 

We propose here a new method to estimate the /3 pa- 
rameter that depends just on the observed monopole and 
quadrupole components of the observed correlation func- 
tion and not on their linear theory representation, includ- 
ing its dependence on cosmological parameters. 

To derive our estimator we start in Fourier space to 
generate as a final product a configuration space expres- 
sion. The observational inputs and their further handling 
belong all to configuration space. Our starting point is 
Equation ( 29 1 , that was shown to be a plausible conve- 



nient desc ription on scales characte rized by quasi-linear 
effects in Percival & White (2009). We also mention 
that such a model is a worth starting point, since it is 
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Fi g. 1 . — Normalized quadrupole (dashed line) defined in Equa- 
tion ( |21| for the streaming model. The linear theory correlation 
function is convolved along the line of sight with a Gaussian ve- 
locity distribution. The asymptote (solid line) Equation ( |22[ l is 
reached at very large scales and furthermore the reading ot the 
plateau around 100 /i~^Mpc is complicated by the presence of resid- 
uals of baryonic acoustic oscillations. 



practically always considered in the literature as a bench- 
mark against which more sophisticated models are con- 
fronted with. Furthermore, quantities in configuration 
space have been shown to be described by perturbation 
theory somewhat better than the Fourier space counter- 
parts. 

The observed monopole and quadrupole can be written 
in terms of the linear spectra: 

Po°^^(fc) = P{k) IJ ^dfik Voifik) (1 + Pf^D' f{<TM 



and 



= /io(/3,a,fc)F(fc) 
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(33) 



==/l2(/3,CT^fc)F(fc) 

The functions h are 

ho{f3, cr^k) = U + ^ + y ) foo{k<Tv) 



(34) 



(35) 



('4/3 4/32 \ 8/?2 
+ ( Y + — I Jo2(k(Tv) + -T^foiikcFv) 



and 



h2W,a,k)^[l+'^ + ^^ f2o{ka,) (36) 
('4/3 4^2 X 

+ ( Y + ^ j /22(fcfT,) + ^/24(fcCT.). 



The functions / are given by 



1 



foiikcTv) = - J d^k fi(^vkfJ,k)'Po{fik)'Pi{^ik) (37) 



and 



f2l{kcFy) = 



where I 



2 

0,2,4. 



d^ik f{crvknk)'P2{^ik)Vi{^ik), (38) 



These functions can be expressed ana- 
lytically in terms of the error or gamma functions for a 
Gaussian distribution or in terms of the arctan function 



for an exponential distribution ( Cole et al.^1995 Peacock 
& Dodds 1996). 



it was shown in Cole et al. (1995) that, in the regime 
where Equation (l29| applies, the ratio 



(39) 



f2°^^(fc) ^ h2{P,Ovk) 
P§^^{k) ho{(3,a,k) 



does not depend on the linear power spectrum P{k) and 
therefore can be used on large scales to determine the /? 
parameter without the need to marginalize on the cos- 
mological parameters. We point out here that the ra- 
tio in Equation ( 39 ) induces also the disappearance of 



any additional function of wavenumber k multiplying 
the power spectrum, possibly approximating non linear 
corrections and/or a scale dependent galaxy bias. Even 
better, the method does not require the knowledge of 
functions of just k in the P°^^{k, iJk) expression, and 
it is therefore more general than the streaming model, 
where P°^^{k, /ife) has to be fully specified. In this work 
we present a new way to exploit these nice properties in 
configuration space, where there might be more favorable 
convergence properties in perturbation theory respect to 
Fourier space. The identities 

V{k) = ho{l3,a,k)P°'°\k)-h2{[3,a,k)P°'^\k) = (40) 

and 

1 1"°° 

D{s) = -— / dkk^j2iks)hoiP,a,k)Pf'ik) 
dkk^j2{ks) h2{f3,a,k)P°''\k) 



1 

iD2{s)~Do{s) 



(41) 



conveniently summarize the basis for our method. We 
first manipulate the ratio in Equation ( 39 1 to reduce it 
to the zero quantity in Fourier space T), and then we ap- 
ply the operator — /q dkk^ j2{ks), with the newly 
formed difference D, that is a configuration space func- 
tion, still being equal to zero. The transformed multi- 
poles are labeled by Z3o and D2. We were then able to 
express the integrals in terms of solely ^q''^, ^2'^'', objects 
in configuration space, and and without any depen- 
dence on the theoretical linear correlation functions. The 
dependence on the cosmological parameters is therefore 
concentrated exclusively on /3. The final result is 
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We now explain all the terms appearing in the differ- 
ence D, that is in general a function of separation dis- 
tance, bu t is equal to zero when the model given by 
Equation ( 29 1 , even thought as multiplied by a generic 
function of k, is valid, the true parameters /3 and 
are used and in the absence of noise. We expect that a 
strong dependence on scale will show up at small scales 
due to highly non linear behavior. T he s ingle and dou- 
ble overline were defined in Equations ( 18 1 and ( 19 ). The 
constants b are 

3 5 



6o=(l 

b2 = 
64- 



4/3 4/32 

8^ 
35 ' 



(43) 



The convolutions C depend on the distribution f{v) and 
on its variance ay. In configuration space, when applied 
on a generic function .g(s), they are 



Coi [9] (s)^ 
C21 [g] (s)^ 

where we used 



diiP2{l-i) 



dvf{v)Vi{fi')gi.s'), 

(44) 



s' = \/s'^ + v'^ — 2fisv 

/IS — V 



(45) 



and Z = 0, 2, 4. The terms with the convolution operators 
C in Equation ( 42 1 are given by 
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^obs 



7=obs' 

5^^ 



(s)= (48) 



1 

2^ 



dkk'j2{k.s)P°^%k)foi{kay) 



and 



C20 [Co°'1 



1 

'2^ 



C22 
1 

"2^ 



(^obs 
SO 



dke n{ks)Pt''{k)f2o{kay) 



iT] is) 



dkk'j2iks)Prik)f22{kay) 



(49) 



(50) 
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(^obs 

so 



5 -^obs 

2?o 



7=ob£ 

2^0 



(51) 



1 

'2^ 



dkk" j2{ks)P°''\k)f2i{k<Jy) 



The formulae in Equations ([23|)-(25) were used to obtain 
the above identities. Furthermore, m Equations (42 1 and 
( 46 ) we employed the definition 



lt\r) EE 
together with the identity 
3a{kr) + j2{kr) = 3 



00 cohs( /\ 



(52) 



°° , , j2{kr') ~ 

dr — = ]2(kr), (53) 



while in Equations ( 42 ) and ( 48 ) we utilized Equa- 
tion ([T9|, 



=obs 5 

^2 = -1^ 



^f^(r')r'4dr', 



(54) 



in conjunction with 

j2(kr) +j4{kr) = 



7= 



j2{kr'y^dr' = -j2(fcO-(55) 



In practical terms it is easy to build a n ef ficient routine 
to compute the estimator in Equation ( 42 1 . The depen- 
dence on P is analytical and only a one-dimensional grid 
on is necessary. A single case with the two parameters 
fixed took us a few seconds on a desktop workstation. 

When ideally applied on the obs erv ed correlation func- 
tions from the model of Equation ( 29 ) , linear theory con- 
volved with a velocity dispersion, the difference estimator 
D completely undoes the convolution with the velocity 
dispersion and returns a zero response if using the true 
P and Gy. This is expected sin ce D was explicitly based 
on the model in Equation ( 29 ) . Much more interesting 
and realistic is when we consider the correlation functions 
built from the Sato and Matsubara iV-body simulations. 
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Fig . 2. — E ffec t of the difference estimator D given by Equa- 
tions BTJ and ||42| l as seen tfirough the two transformed multi polc s 
terms TTq and JJ2 subtracting each other in the rhs of Equation l |41| . 
For better visibility, we have artificially shifted one of the curves 
by 0.7/i~^Mpc and we have multiplied them by the square of the 
separation distance. The correlation functions derived from the 
Sato and Matsubara A''-body simulations described in the text are 
convolved along the line of sight with a Gaussian velocity distri- 
bution according to our prescription and we present an average of 
the different realizations. Since the two curves are very close to 
each other, the prescription D is capable of undoing most of small 
non linearities for scales larger than about 30/i~^Mpc. In this ex- 
ample we fixed the parameters to /3 = 0.85 and = 5.7 Mpc. 
The bulk of the baryonic acoustic oscillations is very similar in the 
two cases. The disparity in the uncertainties, here obtained from 
the diagonal elements of the covariance matrix, comes about be- 
cause the term derived from the monopole has much smaller errors 
respect to the quadrupole one. The role of (T„ consists in trying 
to cope with mild nonlinearities given the range where the fit is 
performed. In the full statistical analysis (t„ is marginalized over. 

It is important to test if and at which scales the action 
of D is able to nullify or markedly reduce the overall ef- 
fect of non linearities and convolution with the velocity 
dispersion. An example of the effectiveness of the pre- 
scription D from Equations (41 ) and (42) is provided in 
Figure[2] It can be noticed that the two terms of the 



right-hand side (rhs) of Equation (41 ), the transformed 
multipoles Dq and D2, are very similar to each other. 
Also, the bulk of the baryonic feature is similar for both 
cases, apart from a small non linear contribution. This 
means that our proposed processing on the monopole and 
quadrupole is capable of bringing the final products to a 
close agreement. The plot is just an example to display 
the general idea of the method. In the next section, we 
perform a full analysis on the simulation data formed 
from the total difference D to find if the best-fitting j3 is 
close to the fiducial value. As already mentioned, we con- 
sider only a Gaussian velocity dispersion. But if needed, 
more effective velocity dispersions might be easily imple- 
mented within the same framework in future analyses. 

Two nice properties of the statistics D are that it does 
not depend on the mean density of the galaxies and that 
it is defined in terms of convolutions. The first attribute 



can be easily proved by looking at Equation (41 ), where 



the presence of the Bessel function assures that the zero- 
mode contribution vanishes. The second feature poten- 
tially provides a certain degree of stability, compared 
with methods that make use of deconvolutions. 

To widen the range of applicability, the method will 
have to be further tested with TV-body simulations that 



include also galaxy bias. However, due to its simplicity 
and practicality, we believe that, with in case necessary 
phenomenological modifications to account for some of 
the effects of non linearities and galaxy bias, our method 
could be a useful tool to study RSDs. 

4. RESULTS 

In a practical situation where data from a galaxy sur- 
vey are available, the quantity D in Equation (42 ) can be 



used directly in configuration space. Here the c orrelation 
functions can be measured, for example, by the |Landy fc| 
Szalay ( 1993 1 estimator. We recommend the application 
of the following function in order to select from the 
data the best /3 and cry. The covariance matrix CovD.y 
relative to the distances Si and Sj for D could be com- 
puted from for example a set of iV-body simulations, like 
we are doing in this work, according to 



1 



irun— 1 

X (i^--(s„/3«^af ) - 5(s„/3*i^fTr)) , (56) 

where irun is the realization index, that reaches the total 
N^nn — 30 in our case. The overbar indicates the mean 
over the realizations and the parameters f3 and (jy are 
considered fixed to fiducial values, that need to be rela- 
tively close to the best-fit values. A few extra iterations 
can update these values if necessary. It is then natural 
to evaluate 



X 



E 

1,3 



D{s,,l3,ay) CoYD-,^D{s,,/3,ay), (57) 



where we have explicitly indicated that Z? is a function 
of distance separation, /3 and ay . We apply this statistic 
to the simulations results. 

The monopole and quadrupole functions furnished to 
us by Sato and Matsubara are evaluated at discrete 
points separated by Ar = 2/i~^Mpc. It turns out that 
the derived covariance matrix CovD^ is highly corre- 
lated and close to singular and in order to be able to 
stably invert it we have used singular value decomposi- 
tion (SVD). The smaller eigenvalues, responsible for the 
instability, are more affected by noise than the larger 
ones and it is necessary to operate a cut upon them, i.e. 
replace their inverses by zeros in the SVD representation 
of the inverse covariance matrix. The situation is more 
conservative for the first eigenvalues but the statistical 
errors are larger, therefore a balance has to be struck. 

We adopt the following method for the eigenvalue se- 
lection. Assuming that the first p eigenvalues surviv- 
ing the cut are almost clean from noise, after having 
found the best-fit /3 and ay, the corresponding '^iH 
have effectively p — 2 degrees of freedom, in the pres- 
ence of a noise close to Gaussian. Hence we choose to 
place the cut where Xp-2 can no longer be accepted as 
a null hypothesis. In practical terms, we summed the 
Xp-2 fo'^ each of the A^run = 30 realizations to form a 



total Xn 
able 



,(p-2) 



statistic. Then we considered the vari- 



2%^ (p-2)' ^^^^ empirically known to approx- 
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imately follow a Gaussian with mean ^y2N^ua{p — 2) — 1 
and unit variance if N^unip — 2) is greater than 30, as it 
is practically always the case in our computations. We 
verified that a stable criterion of rejection is when the 
statistical significance of the test of the correctness of 
the null hypothesis was smaller than 5%. To decide over 
borderline cases we added a Kolmogorov-Smirnov test 
with the same rejection barrier of 5%. Rejection was 
conservatively enacted even if only one of the two tests 
was passed. We found that the two tests gave very simi- 
lar results and that the estimated /?, found by averaging 
over the iVrun =^30 best-fit values obtained by minimiz- 



ing Equation (57 1 for each simulation realization, and its 



error, found from the standard deviation of the best fits, 
were very stable for eigenvalues neighboring the cut. 

Concerning the range of separation distances used in 
the analysis, we fixed the maximum radius to Smax = 
80/i"iMpc, as with our box volume no statistical im- 
provements were evident for larger separations. Instead, 
we varied the minimal separation, that should crucially 
indicate to us where systematics due to non linearity 
are important. The best-fit maximum likelihood val- 
ues of the parameters constitute marginalized samples. 
As shown in Figure[3] we found that for a minimum 
separation such that Smin > 30/i~^Mpc the retrieved 
f3 was compatible at la with the fiducial input value 
f]^'^ = 0.85 with a marginalized statistical error of 5%. 
While for separations such that s,„in > 35/i~^MpCj even 
the errors on the mean were consistent with (3 . At 
Smin = 35 ft.~"'^Mpc wc found that the marginalized statis- 
tical error on f3 is 6%. We conclude that for dark matter 

boxes of volume Vg = I {h~^ Gpc)^, we obtain reliable 
results for Smin > 30/i^^Mpc, since the statistical errors 
are large enough to provide for possible systematic ones. 
While for Smin > 35/i~^Mpc we do not see evidence for 
bias and we could imagine to apply our method to larger 
volumes with smaller statistical uncertainties. 

In order to theoretically estimate the magnitude of the 
error on the parameter /3 given a galaxy survey config- 
uration, we employ the Fisher matrix formalism applied 
to Fourier space. This technique allows us to propagate 
the uncertainties on measured quantities into an uncer- 
tainty on /3. We pursue the following schema: we in- 
terpret the difference I? as a function depending on the 
monopole and quadrupole projections of the observed 
power spectrum, as in Equation ( 40 ) ; we compute the 



errors on those projections; we propagate those errors 
on P; finally we transfer the uncertainties on the /3 and 
Ct, paramet ers. This calculation is based on the model 
of Equation! 29 ) , linear theory convolved with a velocity 
dispersion. We choose the same cosmology used in the 
simulations. 

The variance o f the measurement err or on the observed 
power spectrum ( Feldman et al.|[l994 ) is 



(27r) 



1 



(58) 

where fc is a specific wavenumber at the center of the 
interval Afc, is the survey geometric volume, and n is 
the mean number density of the tracer. We can then 
build the covariance matri x for the multipoles of the 
observed power spectrum (Yamamoto et al. 2006 Ya- 
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Fig. 3. — Estimations of the f} parameter from the Sato and 
Matsubara Af-body simulations. We show the marginaUzed Icr er- 
ror bars for ranges of distance separations going from a variable 
*min to a fixed Smax = SO h~^}Apc, after which there is no statis- 
tical gain, given the simulations box volume. The data have been 
obtained by averaging and forming the standard deviations of the 
marginalized samples produced by a maximum likelihood analysis 
applied to every realization. The fiducial value depicted by a solid 
line is compatible with the data from Smin = SO/i-^Mpc. While 
the data are consistent with no bias from Sniin = 35 h~^Mpc, since 
the errors on the mean, shown as the smaller superimposed error 
bars, start to be less that 2cr from the fiducial value. 




mamoto||2003[ ): 

(2/ + l)(2m+l) 



Gim{k) — ■ 



(59) 



It is possible now to calculate the diagonal covariance 
matrix for the difference operator V of Equation ( 40 ) : 



CovV{k) ^ ^ 



(m{k) dV{k) 

-Ghn{k) (60) 



dP°^^(k) ' 'dP^^ik)' 



The needed partial derivatives are 
&D{k) 



and 



dP§^^{k) 
dV{k) 



ho{l3,ayk). 



(61) 



(62) 



dP^^^k) 

Finally, the Fisher matrix for /3 and is specified by 

''- = E^'=ovi,(.)-?gM. (e) 



dpa 



dpb 



where pa = /?, (Jv indicates the parameter vector and the 
partial derivatives can be taken either numerically or in 
some cases analytically (if for example the distribution 
f{v) is either a Gaussian or an exponential). We use as 
limits of the sum /cmin ~ 27r/(T4)''^^'^'' and utilize a vari- 
able maximum wavenumber k^ax to test the magnitude 
and trend of the errors. 

In Figure[4] we show the percent fractional error for 
marginalizea errors on /3 for a survey centered at red- 
shift z = 1, of volume Vs = 1 (/i~^Gpc)^, full sky. 



9 



«5. 

4- 







25 



30 35 
Smin [Mpc/h] 



40 



Fig. 4. — We show the agreement between the percent frac- 
tional marginalized errors on /3 resulting from the A^-body analysis 
(dashed thick curve) and a Fisher matrix study (solid thin curve) 
carried out in Fourier space and having a maximum wavenumber 
scaled as fcmax = l-Svr/smin- The Fi she r matrix method was based 
on the quantity D from Equation BOjl and refers to the ACDM 
cosmology used in the simulations and is based on a dark matter 
survey at redshift 2 = 1, of volume 1 (/i^^ Gpc)'^, full sky, with 
negligible shot noise. 



with enough dark matter particles to dominate shot 
noise. Our Fisher matrix study is based on Fourier 
space and to relate it to configuration space we need to 
link wavenumbers to distances. We find that the scaling 
fcmax = Svr/j'min givcs thc Fishcr errors of similar magni- 
tude and trend to the errors estimated from the simula- 
tions. While a scaling A:i„ax = l-STr/rmin assures that the 
magnitude and the trends of both the errors are closely 



matched. We note that a similar scaling fc. 



TT/rn 



was chosen by tunin g in Okumura & Jing (2011) and 
Reid & White (2011) to provide a precise link between 
I'burier and configuration space. 

The constraint on l3 can be used to extract information 
on the growth function if the galaxy linear bias is mea- 
sured independently. For instance, the amplitude of the 
matter fluctuations quantified by cr™(2 = 0) and mea- 
sured from a cosmic microwave background radiation ex- 
periment could be first scaled to z = 1, and then the bias 
could be estimated by & = cr|(z = l)/a^{z = 1); the 
superscripts m and g indicating, as usual, matter and 
galaxies. This procedure wa s followed recently, for ex- 



ample, in Guzzo et al. ( 2008 ) , where the fluctuation am- 



plitude derived by the Wilkinson Microwave Anisotropy 



Probe experiment ( Spergel et al.||2007| ) was utilized. 



5. CONCLUSION 

We based our analysis on basic RSD streaming models, 
frequently used as benchmark models in the literature, 
and often shown, for some galaxy groups, to work rea- 
sonably well on quasi-linear scales. 

We have introduced a new method to determine the 
RSD /3 parameter in configuration space on quasi-linear 
scales. The statistics D of Equation (42 ), presents the 
following advantages: it is not necessary to specify the 
linear theory predictions for the correlation functions; 
there is no dependence on the mean galaxy density; only 
convolutions are involved, assuring stability; it does not 
rely on galaxy linear bias, apart from the implicit depen- 
dence within /?. Random errors on the determinations of 
the galaxies redshifts are effectively incorporated in the 
model. 

The method can actually be considered more general 
than the streaming model because it is insensitive to 
a factor multiplying the power spectrum that depends 
only on the wavenumber, which could arise from a quasi- 
linear correction and/or a scale-dependent galaxy bias. 
Furthermore perturbation theory might have better con- 
vergence properties in configuration space, wher e our es- 
timat or operates, compared to Fourier space ([Sanchez 
et al.||2008| [Matsubara||2008[ [Taruya et al.||2009| r 

In addition, we have noticed a baryonic feature at 
about 100 Mpc in the normalized quadrupole in con- 
figuration space, that should not have a major impact on 
the extraction of f3 with our method. 

We have also car ried out analyses based on the TV- 
body simulations of Sato & Matsubara (2011) and the 
Fisher matrix method, finding that errors of a few per- 
cent on /3 are feasible with a full sky, 1 {h~^ Gpc)"^ sur- 
vey centered at a redshift of unity and with negligible 
shot noise and that for minimum separations such that 
Smin > 35/i~^Mpc there is no evidence for a bias in the 
estimation of l3. 

It is going to be crucial to be able to measure and in- 
terpret in detail the growth of fluctuations, because the 
information contained therein could indicate deviations 
from general relativity and constrain the cosmic expan- 
sion history. Since subtleties in the systematics of cosmo- 
logical observations do not guarantee that the theoreti- 
cal equivalence of Fourier and conflguration space can be 
easily reproduced in practical observations, it should be 
convenient to have the widest possible arsenal of meth- 
ods, involving both points of view. 

We are thankful to Joe Silk for reading and com- 
menting the manuscript. We are particularly grateful to 
Masanori Sato for having provided us with the iV-body 
simulations results used in this work and for discussions 
on their features. 
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